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Indirect Field Oriented Control (IFOC) is very popular technique in 
industries due to its simple designing and structure, in comparison to the 
direct method since it requires flux and torque estimators. This paper 
presents detail calculation of all values needed for the development of an 
indirect rotor flux oriented induction motor model fed by a current-controlled 
voltage source inverter. With the calculated values for 1 hp 50 Hz 3 phase 
induction motor; an indirect rotor flux oriented model is developed using 
MATLAB/Simulink. The simulation results are outlined and discussed. 
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1. INTRODUCTION 

Induction motors finds their wide use in industries and homemade applications as it possesses 
several advantages. Despite of various advantages, control of an induction motor is difficult in comparison to 
the DC machines because in DC machine, field and armature fluxes are 90° apart, they do not interact with 
each other and hence flux of the DC motor can be controlled by field current and torque by armature current. 
However, in case of induction motor, stator and rotor fluxes interact with each other due to non orthogonal 
and varying orientations. Hence independent control of flux and torque is difficult. However the difficulty in 
control of Induction motor has been removed by the introduction of several control techniques viz. Scalar 
Control, Vector or Field Oriented Control, Direct torque control etc.[l]-[6]. The Vector or Field Oriented 
Control (FOC) was introduced by Blaschke in 1972. With this strategy, IMs can be controlled like a DC 
machine by maintaining a fixed orthogonal orientation between two fields. 

Field Oriented Control can be implemented in three field or flux orientations as there are three flux 
linkages available in an Induction motor. They are- stator, air gap & rotor flux linkage. Rotor Flux 
Orientation is a technique that provides a method of decoupling the stator currents of squirrel-cage induction 
motor in to two components in the synchronously rotating reference frame (d-q frame):- one, the magnetizing 
current ids*, which is to excite the motor flux and the other torque current iqs*, which is to generate the 
electromagnetic torque respectively. Therefore, it provides independent control of torque and flux, which is 
similar to a separately excited DC machine. Basically, there are two methods of Field Oriented control based 
on the way of determining the rotor angle - Direct method proposed by Blaschke in which rotor angle is 
obtained by the terminal voltages and currents and Indirect method proposed by Hasse in which rotor angle is 
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obtained by rotor position measurement and machine parameters [7],[8]. Rotor flux orientation (RFO) is the 
original and usual choice for the indirect field oriented control . Hence the name Indirect Rotor Flux 
orientation. 

The Field Oriented Control technique needs more calculations than the other control methods. The 
main contribution of this paper is to compute all the values needed for the implementation of indirect rotor 
field oriented control, viz. parameters of the motor, the rated rotor flux, the rated stator d-axis and q-axis rms 
and peak values of currents and the value of the slip speed for rated torque operation, selection of the DC 
voltage etc [9],[10],[11]. It is also necessary to determine the parameters of the speed PI controller-Kp & Ki. 

This paper is organized as follows: Section 2 presents the calculations for parameter determination 
of the Induction motor under test, Section 3 describes the calculations required for design of Indirect Rotor 
Flux orientation control scheme. Section 4 present the SIMULINK model development for Indirect Rotor 
Flux oriented Induction motor using MATLAB. In section 5 simulation results & discussions are presented 
and the paper ends with the conclusion. 


2. CALCULATIONS FOR PARAMETER DETERMINATION 

Determining the parameters of Induction Machines is very important step while considering Field 
Control Technique. The parameters of the induction motor are determined by performing no-load and 
blocked rotor test [12]. The name plate data of the induction motor under test is mentioned in Table 1. 


Table 1 Name plate data of Induction Motor under test 


Nominal 

Power 

(Watts) 

Voltage (line-to- 
line) 

(Volts) 

Frequency 

(Hertz) 

Inertia 

(Kg.m A 2) 

Friction 

factor 

(N.m.s) 

No.of pole 
pairs 

Speed 

(rpm) 

Rated 

Current 

(Amp) 

Pn 

Vn 

fn 

J 

F 

P 

Nn 

Isn 

750 

415 

50 

0.005776 

0.00328 

2 

1490 

2 


The results of the no-load and blocked rotor tests on 1 hp, three-phase, 415 V, four-pole, 50 Hz, 
Squirrel-cage induction motor are given in Table 2 and Table 3 


Table 2 No-load Test data 


Line to line Voltage 
(Volts) 

Line Current 
(Amp) 

Input Power (watts) 

Speed (rpm) 

Frequency (Hertz) 

Vnl 

Ini 

Pnl 

Nnl 

Fnl 

412 

1.4 

144 

1497 

50 


Table 3 Blocked Rotor test data 

Line to line Voltage 
(Volts) 

Line Current 
(Amp) 

Input Power (watts) 

Speed (rpm) 

Frequency (Hertz) 

Vbr 

140 

Ibr 

2.58 

Pbr 

386 

Nbr 

0 

Fnl 

50 


Stator resistance 9.395 O per phase. From No-load test data, following parameters are computed. No load 
Voltage per phase: 


t / u Vnl 

Vph =s 

No load Power per phase: 


( 1 ) 


Pph=™ 


( 2 ) 
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No-load power factor: 


cos (/nl - 


Pph 


Vph * Ini 

Magnetising current: 

Im= Ini* 

Core loss current: 

Ic = Inl*cos(tfnl) 
Magnetising or Mutual Inductance: 


Lm= - 


Vph 


2pi* fnl*lm 


From Blocked rotor test data, following parameters are calculated- 
Short circuit voltage per phase: 


Vbr 

VSCpk- —pr- 

V3 


Short circuit power per phase: 


Pscph=^f 


Short circuit power factor: 


cos (j)SC - 


Pscph 

Vscph*Ibr 


Short circuit impedance: 


Zsc 


Vscph 

Ibr 


( 3 ) 


(4) 

(5) 


( 6 ) 


(7) 


( 8 ) 


(9) 


( 10 ) 


Rotor resistance: 

Rr=Zsc*cosl$sc)-Rs (11) 


Equivalent reactance: 

Xeq= Zsc* smic/jsc) = XI s+ Xlr 1 


( 12 ) 


Where Xls= Stator leakage reactance 

Xlr’= Rotor leakage reactance as referred to stator 

The connection between Xls and Xlr’ for the most common types of industrial Induction motors 
grouped by the National Electrical Manufacturer’s Association (NEMA) is shown in Table 4.[11] 


Simulation of an Indirect Rotor Flux Oriented Induction Motor Drive... (Vivek V. Puranik) 



1696 n 


ISSN: 2088-8694 


Table 4. NEMA Class motor classifications 


NEMA Class 

Connection 

Class A and D 

Xls=Xlr’ 

Class B 

Xls=(2/3)Xlr’ 

Class C 

Xls=(3/7)Xlr’ 


Since most of the general purpose motors are Nema class B, for Nema Class B Induction motors, 

Stator Leakage Reactance! 

Xls=0A*Xeq (13) 

Rotor Leakage Reactance as referred to stator: 

Xlr 1 = 0.6* Xeq (14) 

A .m file is created containing all the data mentioned in Table 1 to Table 3 and all the formulas 
represented by (1) to (14) to calculate the parameters of the Induction motor. Thus the calculated parameters 
of the induction motor under test are tabulated in Table 5. 


Table 5.Parameters of Three Phase Induction Motor 


Stator Resistance 

Stator Leakage 

Rotor Resistance 

Rotor Leakage 

Mutual Inductance 

(Ohm) 

inductance 

(Henry) 

(Ohm) 

inductance 

(Henry) 

(Henry) 

Rs 

Lis 

Rr’ 

Llr’ 

Lm 

9.395 

0.0350 

10.444 

0.0525 

0.5492 


3. PROPOSED METHODOLOGY AND ITS DESIGN 

The vector control technique chosen here is an indirect rotor flux orientation. The basic principle of 
the indirect rotor field oriented control scheme of an induction motor is to lock the rotor flux vector with d-q 
co-ordinate reference frame. These results in decoupling of stator currents which controls the flux and the 
torque separately by controlling stator direct axis current ids and stator quadrature axis current, iqs 
respectively like a separately excited DC machine. Locking of d-q co-ordinate reference frame revolving 
with rotor flux needs the information about the magnitude and the position of the rotor flux. Lig.l shows the 
block diagram of the Indirect Rotor Lield oriented control scheme. The torque producing component of the 
stator current, iqs* is: 


iqs* = 



Te* 

\\f/r\est 


(15) 


where Lr is the rotor inductance, Lm is the mutual inductance, and |\|/r|est is the estimated rotor flux linkage 
given by: 


I TV I est = L ,n id 

1 + TT.S 


where x r = Lr / Rr is the rotor time constant. 

The flux producing component of stator current ids* is: 


ids* = 


I TV I * 
Lm 


(16) 


( 17 ) 
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The rotor flux position 0 required for coordinates transformation is 


6 = \(coe + cosl)dt. 

The slip frequency cosl is: 


U71 Lm Rr . 

Wsl =- xqs 

I Wr I est Lr 


(18) 


(19) 



Figure 1 Block diagram of Indirect Rotor Field Oriented Control Scheme 


The q and d-axis voltages in synchronous reference frame are [3] given by: 

Vn*J 2 

V* = ° 

The induction motor model in electrical synchronous reference frame is given by: 
Vq =(Rs + Lsp)i q +o) L s i d +L m pi qr+(D L m i dr 

v ds =~co s L s i qs + ( Rs+Ls pH ds -co s L m i q , +L '»pi dr 

Vqr = L mPi qs + { COs ~ COr^L + (Rr+ ^ h,r + C CQ S ~ COr^L 

V d =-(CQ~CO) L m i q +L m pl d +(cO-CO) L r i q +(Rr+Lrp)r 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 
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The steady state currents are thus evaluated by using synchronous reference frame equations with the 
substitution of p=0 and the slip speed being zero. The stator and rotor flux linkages in electrical synchronous 
reference frame are given by: 


UJ — J / _L J i 

T qs lj s L qs ' Ll m L qr 

(26) 

ds Lsids ^m^dr 

(27) 

Ipqr ~ L r iqr T L m iq S 

(28) 

Ip dr ~ L r i(ir 3" L m ids 

(29) 


Substituting the corresponding values, the resultant stator & rotor fluxes are given by: 

'¥s = il%+'¥l 

Crated 

3.1 Calculation of Rated Torque 

Under no-load & rated condition, rated torque is given by: 

T L =Te=Tsh= Pn/co n ~5 N-m. 


(30) 

(31) 


(32) 


All the above formulas from (15) to (32) are included in the earlier created .m file and the commanded values 
of voltages, currents and fluxes are computed. 

3.2 Speed Controller 

The most commonly used controller for the speed control of induction motors is conventional PI controller. 
The transfer function of PI controller is given by: 


G c (s) = 


K p s + K t 
s 


(33) 


Proportional gain and the integral gain, Kp & Ki, are the parameters of PI controller. These parameters are 
determined using manual tuning technique. In manual tuning process, initially Kp was kept equal to 1 & Ki 
equal to 0.0001. Then, Kp was increased till the output shows some oscillations. Then Ki was changed to 
remove steady state error. 

3.3 Selection of Dc Voltage (V DC ): 

The induction motor is fed by a current-controlled PWM inverter which has 2- voltage levels. For a 2-level 
voltage, line-to-line voltage is: V L _ L =0.612 *m index *Vdc 
Since modulation index < 1, DC link voltage is: 


^DC - 


Yl=l 

0.612 


So the inverter is built with Vdc= 700Volts. 


(34) 


4. DEVELOPMENT OF SIMULINK MODEL 

Based on the equations presented in the section 3 and using various tool boxes available in the 
SIMULINK library, model is developed for each block. Step by step development of model is presented here. 
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4.1 Speed Regulator Block 

As mentioned in section 3.2, the commonly used speed regulator is conventional PI controller. The simulink 
model of it is implemented as shown Figure.2. 



Figure 2 Simulink model of Speed PI controller 


4.2 IFOC Block 

This block is implemented as shown in Figure 3. It consists of 7 main sub-blocks: Ids* computation, Iqs* 



Figure 3 Simulink model of IFOC 


Figure 4 Simulink model of Theta 
Computation 


computation, ABC-to-dq transformation, dq-to-ABC transformation, Flux computation, Slip computation and 
Theta computation blocks. Simulink model of some of these blocks are presented in the following subsection. 
Based on (18), Simulink model of this block is implemented as shown in Figure 4. 

4.2.1 ABC-TO-DQ Transformation 

This transformation is called as Park’s transformation. A Simulink model is developed based on the 

(35). 


V 

_ 2 

cos(0) cos(0-;k) cos(# + x) 

4 

J 

_v 

~~ 3 

-sin(0) -sm(0Fy) 

L b 

Jc_ 
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CZ> 

Id 

G> 

lb 


d> 

lc 


ci> 

Theta 


(cos(u(4 -)T u(1 )+■ cos(u(4>2* p \S3J u(2 )+■ cos(u(4)4-2 K p i/3 T u(3)) 


-(si nCuC4))T u(1 )*si n(u(4>2* p if3? u(2)f si n(u(4)4-2* p i/3)f u(3)) 


F cnl 


Odin 


Figure 6 Simulink model of ABC-to-dq transformation 


4.2.2 DQ-TO-ABC Transformation 

This transformation is called as Inverse Clarke’s transformation. A Simulink model is developed 
based on (36). 



Figure 7 Simulink model of dq-to- ABC transformation 


4.3 Current Regulation 

Current regulator used in this field oriented control is hysteresis band current controller. It is one of 
the several types of PWM techniques. It is basically a feedback current control method of PWM where the 
actual currents are compared with the reference currents within the specified hysteresis band. Thus the inputs 
to the Hysteresis controller are three phase current errors and the outputs are the switching pulses to the 
inverter. Fig. 8 below shows the Simulink model of Hysteresis current controller. Simulink model for other 
blocks viz. Iqs* Computation, Flux computation and Ids* computation are implemented as per (15) to (36). 
Fig.9 shows the complete simulink model of the indirect rotor field oriented Induction motor which is now 
ready for simulation. 


5. SIMULATION RESULTS AND DISCUSSION 

Before starting the simulation, a .m file prepared earlier is executed at the Matlab prompt so that all 
the values of the variables in the Simulink model will be assigned to them. A switching frequency of inverter 
is taken as 20KHz and a DC link voltage is taken as 700 V. 
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Figure 8 Hysteresis Band current regulator 





Figure 9 Simulink model of Indirect Rotor Field Oriented Induction Motor 


Initially the reference speed is set to a constant value of 100 rad/s (955rpm) and the load torque is 
also kept constant at 0 N-m. Setting the correct parameters of the speed PI controller and hysteresis band is 
an important task in simulating the model for getting the accurate results. The values of Kp, Ki and hysteresis 
band are to be adjusted number of times and simulation is run every time until accurate results are 
obtained. There are many frequency domain and time domain methods available for designing of a PI 
controller including the most popular Ziegler-Nichols method of tuning PID controller. But all these methods 
do not guarantee that the controller parameters found by them gives you accurate performance 
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characteristics. They just give the initial guess for the controller parameters and finally you have to tune the 
controller manually. 

Here the speed PI controller is tuned manually. The simulation is run for the 3 sec and the responses 
are observed on the respective scopes. The response curves of Rotor speed, electromagnetic torque (Te*), 
flux, stator current (Iabc) vs time are observed and are shown in Figure 10 to Figurel3. Figures on the left 
hand side are for Kp=l , Ki=0.001 and h=l whereas those on right hand side are for Kp=4 , Ki=0.15 and 
h=0.006. After comparing the figures on the left hand side with right hand side, it is cleared that due to 
improper selection of PI controller parameters, speed is not reaching to the set speed and hence the flux is 
also not reaching to the calculated value of 1.012 Wb. Moreover due to improper selection of hysteresis band 
width h, disturbance in the torque and stator currents are clearly observed. Simulations are carried out for 
various reference speeds. 




(b) 


Figure 10 Rotor Speed vs Time 



(a) (b) 

Figure 11 Electromagnetic Torque vs time 
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(b) 


Figure 12 Rotor Flux vs time 



(a) (b) 

Figure 13 Stator currents Iabc vs time 


6. CONCLUSION 

Simulink is a graphical extension to MATLAB to develop a model and then perform experiments on 
it for predicting the behaviour of the systems over time without building it practically. So it is a first step 
towards designing a system. In this paper detail calculations necessary for the development of the Indirect 
rotor field oriented control of induction motor drive using Matlab/ Simulink has been presented .The model 
so developed can be used to study the dynamic behaviour of the induction motor by applying step change in 
the reference speed and the load torque. From the present model developed, one can implement various types 
of control schemes by adding subsystems in the present model. This paper will be helpful for the post¬ 
graduates and researchers in real time implementation of the drive. 
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